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Abstract 

We present a study of dissociative electron attachment and vibrational excitation processes in 
electron collisions with the CF3CI molecule. The calculations are based on the two-dimensional 
nuclear dynamics including the C-Cl symmetric stretch coordinate and the CF3 symmetric de- 
formation (umbrella) coordinate. The complex potential energy surfaces are calculated using the 
ah initio K-matrix method. The results for dissociative attachment and vibrational excitation of 
the umbrella mode agree quite well with experiment while the cross section for excitation of the 
C-Cl symmetric stretch vibrations is about a factor of three low as compared to experimental data. 
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12 I. INTRODUCTION 



13 Dissociative electron attachment (DEA) to polyatomic molecules typically involves mul- 

14 tidimensional nuclear dynamics. However, because of big computational work necessary 

15 to obtain multidimensional complex {i.e., including both real and imaginary parts) energy 

16 surfaces, most of theoretical DEA calculations were performed in one-dimensional approx- 

17 imation. In these calculations it is usually assumed that the DEA process involves one 

18 reaction (dissociating) coordinate, roughly corresponding to one of the normal modes of the 

19 target molecule. This approximation is sometimes too crude, and sometimes completely un- 

20 justified. Therefore, a lot of effort was devoted recently to calculations of multidimensional 



21 DEA dynamics 



y-[6|. These calculations address two important problems in the physics of 

22 DEA processes. First, we want to know which dissociation channels are the most important 

23 and what is the energy range where a particular bond breaking can occur. This informa- 

24 tion is especially important for chemical control. Second, we want to know the importance 

25 of different vibrational modes in a particular DEA process and the final-state vibrational 

26 energy distribution in the fragments. 

27 The most common method for studies of dynamics on multidimensional surfaces is the 

28 wave-packet propagation technique 7|. Recently, this approach was used to calculate the 

_or HoO j^l, employing 



10 



An alternative 



29 DEA cross sections for several polyatomic molecules, e.g. for CO2 ||8i_ 

30 the multi-configuration time- dependent hartree (MCTDH) method 

31 is to develop classical and semiclassical methods for treatment of nuclear dynamics. We 

32 recently reformulated ^ the quantum method for DEA in terms of the time-independent 

33 Schrodinger equation and connected this treatment with the classical approximation. The 

34 method was applied to the process 

e + CF3a(z/2, z/3) CFsCr ^ CF3(z/^) + CI". (1) 

35 Here z/2 and stand for the symmetric deformation vibrations (so-called "umbrella" mode) 

36 and the symmetric stretch vibrations, z/g represents the umbrella mode of the free CF3 radi- 

37 cal. The two-mode approximation for this process can be justified by existing experimental 

38 data [12] on vibrational excitation (VE) of this molecule. The two-dimensional potential 

39 energy surface (PES) was calculated ab initio. However, we used a model semiempirical 

40 width (the imaginary part of the complex PES). This led to some inconsistencies and in- 

41 stabilities in our calculations discussed below. In the present paper we employ the ab initio 



42 molecular R-matrix method for calculation of the complex PES. This allows us to remove 

43 the deficiencies in our previous calculations and improve agreement with experimental data. 

44 We also calculate VE cross sections for the C-Cl symmetric stretch and the umbrella mode. 

45 The rest of the paper is organized as follows. In Sec. [Ill we discuss construction of the 

46 complex PES from the ab initio calculations and the our theoretical approach to treatment 

47 of nuclear motion. In Sec. IIIII we discuss details of our R-matrix scattering calculations in 

48 the fixed-nuclei approximation. In Sec. IV we present our results for DEA and VE. 



49 II. THEORETICAL APPROACH 

50 In our previous work [ll] we constructed a two-dimensional local complex potential (LCP) 

51 model for DEA and resonant VE of CF3CI. The model takes into account the C-Cl stretching 

52 vibrational mode and umbrella vibrational mode of the CF3 radical. The LCP model [l| was 

53 based on a complex PES constructed from the one-dimensional potential curves using an 

54 arbitrary extension in the coordinate corresponding to the umbrella vibrations of the CF3 

55 fragment. In addition, the real part of the complex potential curve for the temporal anionic 

56 complex was obtained using the ab initio methods for bound-states calculations, while the 



57 imaginary part was obtained by fitting 



13| the experimental results of Mann and hinder 



58 [12]. Cross sections of the DEA calculated using the two-dimensional model are factor of 

59 three higher than experimental values and one-dimensional non-local calculations [l^. This 

60 discrepancy was attributed to the inconsistencies between the real and imaginary part of 



61 the complex PES used in the model [l]. It is the aim of this work to construct a two- 

62 dimensional model of the nuclear dynamics along the same lines as in pj], however using 

63 a more consistent complex PES. We performed molecular ab initio R-matrix scattering 

64 calculations in the fixed-nuclei approximation for a set of nuclear geometries including both 

65 degrees of freedom corresponding to the C-Cl distance R and F-C-Cl angle {}. We fitted 

66 the eigenphases obtained from these R-matrix calculations at energies close to the resonance 

67 position to the Breit-Wigner formula with an energy- dependent background and using 

68 the resonance position and width we constructed the complex PES. This approach is free of 

69 any presumption on the dependence of the complex PES on the F-C-Cl angle and removes 

70 the need for any arbitrary extension in this coordinate used in In addition, the fixed- 

71 nuclei resonance width is calculated at the same level of the theory as the position. 



72 The ab initio molecular R-matrix method is well known and widely used for fixed-nuclei 

73 calculations of electron collisions with small and medium-sized molecules. We refer the 



74 reader to a recent review article by Tennyson 15| and to references therein for description 

75 of the method and its technical implementation. 

76 The eigenphase sums calculated using the R-matrix method for a set of nuclear geome- 

77 tries were fitted using the Breit-Wigner formula taking into account the dominant dipole 

78 component of the potential in the outer region. The Breit-Wigner formula is equivalent 

79 to the one-pole approximation to the R-matrix [l^. This approximation is based on the 

80 assumption that at the energies close enough to the resonance position the R-matrix can be 

81 well approximated by the following expression: 

RiE) = R, + (2) 

82 where Rq is the background R-matrix including all the terms due to remaining poles of the 
33 R-matrix, Ex is position of the pole closest to the resonance, 7a is corresponding amplitude 
84 and E is the scattering energy. To obtain the Breit-Wigner formuula from Eq. ([2]), we 

35 assume that the term Rq is a slowly varying function of energy as it is sum over all the other 

36 R-matrix poles. The eigenphase sum 6 may then be expressed [l^ as 

m = tan-i L ^ ) - (3) 



Ea + Aa-^, 

87 where Fa is the resonance width and Aa is the level shift (amount by which the resonance 
33 energy is shifted from the pole E\). The first term of Eq. ([3]) describes the resonance contri- 
39 bution and the second {(f){E)) the potential scattering contribution. The relations between 

90 Fa, Aa, and general solutions of the Schrodinger equation on the R- matrix boundary are 

91 given in Ref. [l4|. In order to fit the width and the position of the resonance using this model 

92 we first solve the Schrodinger equation with the dipole potential in the outer region. The 

93 corresponding dipole moment for every nuclear geometry is obtained from the ab initio cal- 

94 culation of target properties as a part of the fixed-nuclei scattering calculations. Having the 

95 solutions on the R-matrix boundary we can establish the relation between the background 

96 phase shift (f){E), width Fa, level shift Aa and quantities in Eq. This allows us to fit 

97 the model R-matrix amplitude 7a, pole Ex and constant background R-matrix Rq using the 
93 non-linear least-squares technique to the ab initio eigenphase sum. These directly determine 
99 the width Fa and resonance energy Ex + Aa. This fitting allows for a construction of the 



100 complex PES U{R, r) — ir{R, r)/2 in the region where the anionic state is metastable. In the 

101 region where the resonance turns into a bound state, the corresponding bound state energy 

102 was calculated as well as the potential energy curve V{oo,r) of the free CF3 fragment. The 

103 complex PES was constructed using cubic splines in two dimensions. In order to study the 

104 final-state interaction on the anionic surface during the DEA process, it was necessary to 

105 have the bound anionic PES also for large C-Cl internuclear separations, where the ab initio 

106 results are not available. It was constructed by the extrapolation of the ab initio results 



107 as described in 



l| to satisfy the condition U{R — )■ oo,r) — > V{oo,r). The bound part of 



108 the anionic potential energy surface was extrapolated to match the CF3 fragment potential 



109 curve asymptotically, as discussed in 

no The LCP calculations of the DEA and VE presented here were performed in a similar 

111 way as described in Ref. pj]. The basic equation of the LCP theory reads 

[T, + Tr + U{p, r) - ir(p, r)/2 - E]xs(p, r) = Vdk{p, r)Up, r), (4) 

112 where p and r are the reaction coordinates introduced to decouple the two-dimensional 

113 operator of the nuclear kinetic energy U{p,r) — ir(p, r)/2 is the complex PES of the 

114 temporal anion, Vdkip,r) = >yT{p, r)/27r is the amplitude for electron capture into the 

115 resonance state and Ci{pi^) is the vibrational wave function of the neutral molecule in the 

116 initial state. Tp + % is the operator of nuclear kinetic ener^ corresponding to our two- 

117 dimensional model as discussed in [li]. In our previous work [l] the wave function XE{p,r) 

118 was expanded in the basis of vibrational states 0;y(r) of the CF3 fragment in the harmonic 

119 approximation given by equation 

[T,, + y'^(oo,r)-e,]</),(r) = 0, (5) 

120 where e,^ are the corresponding eigenenergies e^, = De + ul{u + 1/2), ul is the harmonic 

121 frequency of the CF3 radical umbrella mode and V^{oo,r) is the corresponding free CF3 

122 radical potential curve in the harmonic approximation with the minimum corresponding to 

123 the C-Cl bond dissociation energy. As it is explained in projection of Eq. (jlj) on 0^(r) 

124 then yields a set of coupled differential equations for channel wave functions of the variable 

125 p with coupling potential 11^,^1 (p) given by the equation 

U..'{p) = J Mr) [U{p,r) - ir(p,r)/2 - l^'^(oo,r)] M{r)dr. (6) 



126 One note should be made at this point regarding the asymptotic behavior of the couphng 

127 potential. The PES in [1] was extrapolated in such way that U{R -t- oo,r) -t- V{oo,r). 

128 The extrapolation asymptotically matches the ab initio potential energy curve of the CF3 

129 fragment, not its harmonic approximation. As a consequence, limp_^^U^i,r{p) 7^ 0, as can 

130 be seen in Eq. Q. Since the coupling of different vibrational channels of the CF3 radical 

131 doesn't vanish for p — )■ 00, population of different vibrational states of CF3 produced by the 

132 DEA process does not converge well. This deficiency is corrected in the present work. The 

133 harmonic approximation for the potential energy curve V^{oo,r) of the CF3 radical is not 

134 used and the full ab initio potential curve V{oo,r) is employed in both Eq. dSD and Eq. (|H]). 

135 Therefore, in our present calculations we expand the wave function XE{p,r) in the basis set 

136 of eigenfunctions (pu{r) given by the equation 

[Tr + Vioo,r)-e,]ip,ir)=0, (7) 

137 where e,^ are the vibrational energies of the CF3 fragment without the harmonic approxima- 

138 tion. The corresponding coupling potential has the following form: 

Uuu'{p) = j M'r)[Uip,r) -iT{p,r)/2-V{oo,r)]ip,,{r)dr. (8) 

139 Therefore, in the present work the extrapolated PES is consistent with the asymptotic 

140 treatment of the CF3 fragment that allows for convergence of the population of vibrational 

141 states of CF3 produced by the DEA process. 

142 In the present work we also use a different method to solve the system of coupled ra- 

143 dial equations. In Ref. pj] we used direct outward integration of the system of differential 

144 equations from the inner region and inward integration from the asymptotic region with 

145 subsequent matching of the solutions to satisfy the boundary conditions. The direct inte- 

146 gration has a limitation in number of channels included in the calculation. With increasing 

147 number of closed channels included, their exponentially increasing contribution starts to be 

148 more pronounced and the calculation becomes unstable. 

149 Here we employ the multichannel version of the exterior complex scaling (ECS) method in 

150 the discrete variable representation (DVR) basis set [l^, which is free of this problem. This 

151 method has previously been successfully used in the context of nonlocal resonance model 



152 for calculations of DEA and VE in case of diatomic molecules 



17|. 



153 III. SCATTERING CALCULATIONS IN THE FIXED-NUCLEI APPROXIMA- 

154 TION 



155 In order to obtain the two-dimensional complex PES necessary to construct the LCP 

156 model, we performed the R-matrix scattering calculations in the fixed nuclei approximation 

157 for a set of nuclear geometries important for the DEA and VE. For every nuclear geometry we 

158 calculated the eigenphase sum in the energy interval around the fixed-nuclei resonance and 

159 fitted it to the Breit-Wigner formula with background as described above. The resonance 

160 position and width as a function of the nuclear coordinates represent the complex PES 

161 used in the local complex approximation. In the region where the CFsCl^ anion is stable 

162 against autodetachment and the potential becomes real, the surface is represented by the 

163 bound-state energies calculated using the R-matrix approach. 

164 As it was discussed previously in our calculations the PES V and U for the neutral 

165 molecule and for the anion are represented using two coordinates: C-Cl internuclear sepa- 

166 ration R and the distance between the C atom and the plane formed by the fiuorine atoms 

167 r = —Rcpcos^, where ^9 is the F-C-Cl angle and Rcf is the F-C bond length. Since we 

168 do not include the C-F stretching mode into our considerations, Rcf is fixed and set to 

169 the value 1.342 A corresponding to the equilibrium geometry of the neutral CF3CI. The 

170 fixed-nuclei R-matrix calculations were performed for a two-dimensional region of nuclear 

171 coordinates with R from 3ao to 12 ao and with -i? from 55° to 90°. In the present work we 

172 consider excitation of the low vibrational states of the CF3 fragment only, therefore we don't 

173 take into account any effects caused by the flipping of the radical > 90°). 

174 Calculations were performed using Gaussian type orbitals (GTOs) and the UK polyatomic 

175 R-matrix code The highest symmetry available in the polyatomic code is Cg which is 

176 an abelian subgroup of the true Cs^ symmetry of CF3CI. 



177 A. Target representation 

178 The CF3CI was represented using Hartree-Fock (HF) molecular orbitals (MOs). In order 

179 to construct a target model sufficient for the purpose of the dynamical calculations, we 

180 performed several tests with different target models to select the best compromise between 

181 the quality of the target representation and the computational tractability of the (iV -|- 1) - 



182 particle problem. 

183 The CF3CI target states were represented using a complete active space (CAS) config- 

184 uration interaction (CI) wave function. The CF3CI molecule contains 50 electrons. Only 

185 18 of them belong to the inner shells, 32 remaining electrons form the valence shells and 

186 in principle can contribute to the chemical bonds. This complicates the construction of the 

187 target CI model in several aspects. Enough valence electrons should be included in the CI 

188 active space to treat the electron correlation properly. In addition, it is necessary to treat 

189 the target symmetrically and to ensure that all elements of the degenerated pairs of MOs 

190 are included in the CAS. On the other hand, inclusion of each orbital occupied in the HF 

191 ground state leads to a rapid increase in dimension of the target and anionic CI Hamiltoni- 

192 ans. In addition, it is the aim of our fixed-nuclei calculations to study the dependence of the 

193 resonance on the nuclear geometry. This raises further the limitations on the N + 1-particle 

194 CI calculation as we need to repeat it many times. 

195 The first CAS CI target model considered in our calculations includes eight active elec- 

196 trons which occupy four orbitals in the HF ground state determinant with the highest orbital 

197 energies. We allow these electrons to occupy five lowest virtual orbitals (VOs). Using the 

198 notation of the Cg point group, this CAS model can be expressed as follows: 

(la' . . . 15aO'°(la" . . . 6a")''(16a' . . . 21a', 7a" . . . 9a")^ (9) 

199 or 

(la' . . . 15a')=^°(la" . . . 6a"Y\l6ei' . . . 20a', 7a" . . . 10a")^ (10) 

200 as the ordering of VOs changes with nuclear geometry. Using the notation of the C^^ point 

201 group, the active electrons are taken from the following set of HF orbitals: (lci2; lOai, 7e). 

202 The degenerated e-orbitals in this set are well localized on the chlorine atom, the lOai- 

203 orbital is spread along the C-Cl bond and the la2 orbital contributes to the C-F bonds. 



204 The calculations were performed using the 6-311G* GTO basis set |18l. Il9|. 

205 Representation of the target in the subsequent scattering calculation can be qualified by 

206 comparison of several properties calculated using our model with the previously published 

207 results. To check our target representation, we calculated the vertical excitation energy 

208 (as the excited target states are used in the scattering calculations), dipole moment of the 

209 target ground state (as it represents the major contribution to the interaction with the 

210 projectile in the outer region) and static dipole polarizability which allows us to estimate 

8 



211 the representation of the polarization effects in the scattering calculation. Comparison of 

212 values calculated using our models (at equilibrium nuclear geometry) with the previously 

213 published data is shown in Table HI As it can be seen there, this target model (denoted 



TABLE I. Target properties calculated using different CAS CI models at equilibrium nuclear ge- 
ometry and their comparison with the previously published results. 



214 
215 





Model 1 


Model 2 


Bibliographical data 


GTO basis 


6-311G* 


cc-pVDZ 




# of active electrons / ^ of VOs 


8/5 


4/11 




Ground state energy (a.u.) 


-795.787343 


-795.627044 




Vertical excitation energy (eV) 


9.545 


9.057 


7.7 ±0.1 [20, 21j 


Vertical excitation energy '^E (eV) 


8.476 


7.897 




Dipole moment (a.u.) 


0.378 


0.410 


0.197 [22] 


Static dipole polarizability (a.u.) 


1.32 


0.158 


38.6 [22] 


TABLE II. Vibrational frequencies calculated using 


different CAS CI models and their comparison 


with the previously published results. 








Model 1 


Model 2 


Bibliographical data 


UJ2 (cm ^) 


854.42 


956.48 


775.12 [23j 


UJ3 (cm^i) 


433.08 


534.13 


463.33 [23] 


CF3 fragment (102) (cm^-*^) 


755.53 


745 


701 [24j 



216 
217 



218 as Model 1 there and elsewhere in the text) gives a reasonably good representation of the 

219 ground state dipole moment, as it is relatively small, although our calculation gives us a 

220 larger value than found in the experiment. The vertical excitation energy calculated here is 



221 almost 2eV higher than the experimental value [2l|]. It suggests that the representation of 

222 the target excited states is limited. According to our knowledge, there is no experimental or 

223 advanced theoretical calculation of the lowest excited state available. Our target model 

224 gives its energy 8.453 eV above the ground state. In our scattering calculations, we are 



225 interested in scattering energies below 4 eV, where all the electronically excited channels are 

226 closed. We expect (and our test calculations described below suggest it) that inclusion of 

227 the low excited target states in the close-coupling (CC) expansion leads to a small correction 

228 of the resonance position and width only. 

229 In order to estimate how well are the polarization effects represented in our scattering 

230 calculations, we evaluated the static dipole polarizabihty of the target ground state. We 

231 used the sum-over-states formula with the set of target states included in the scattering 

232 calculation. We found that the lowest six target states represent a considerable contribution. 

233 Adding more target states did not increase the value. As can be seen in Table [U our 

234 calculated value is considerably smaller than the experimental results. Since the higher 

235 excited states do not have significant contribution to the polarizabihty, this discrepancy 

236 can be possibly attributed to low number of active electrons and orbitals used in our target 

237 model. Good representation of the polarizabihty would require presence of more diffuse MOs 

238 which are too high in energy to be included in the present target calculations. However, the 

239 previously published one-dimensional calculations of the nuclear dynamics for CF3CI [13,] 

240 show that the polarization effects have a minor effect on the results of the LCP calculations. 

241 Target properties discussed above are all calculated at the equilibrium nuclear geometry. 

242 However, for our calculations of the nuclear dynamics it is important to know how well can 

243 our target model reproduce the harmonic vibrational frequencies of the neutral molecule. 

244 Comparison of our calculated values with the previously published result is showed in Ta- 

245 ble [TTl Our target model 1 gives the values with 10% accuracy when compared with the 

246 experimental results. This small difference can be possibly attributed to the neglect of other 

247 vibrational degrees of freedom in our model. 

248 We can conclude that our CAS CI model 1 represents the neutral target sufficiently 

249 enough to be used in our ab initio K-matrix calculation. R-matrix calculations performed 

250 at the equilibrium nuclear geometry can be used to obtain the vertical attachment energy. 

251 Correct value of this quantity is essential for our calculations of the resonant nuclear dy- 

252 namics. However, the CAS CI model 1 of the target leads to the value which is 1.3 eV 



25|. Additional testing 



253 above the experimental value determined by Aflatooni and Burrow 

254 calculations show that it is due to insufficient treatment of the electron correlation in the 

255 CC expansion of the N + 1-electron wave function. Since further increase of the number of 

256 active electrons or VOs included in the target model leads to intractable diagonalization of 

10 



257 the + 1-electron Hamiltonian, we decided to model the electron correlation by modifica- 

258 tion of the primary GTO basis set, where we modified the gaussian exponent corresponding 

259 to the p-orbital of the chlorine atom in order to get the energy of the lowest unoccupied 

260 molecular orbital (LUMO) closer to the experimental value of the vertical attachment en- 

261 ergy. We expect that this modified LUMO helps to represent the discrete component of the 

262 resonant wave function of the temporal anionic complex better than the linear combination 



263 of VOs obtained in model 1. This idea is similar to the method developed in Ref. 



26| 



264 used to calculate the position and width of the resonance. The modified HF orbitals were 

265 calculated using the cc-pVDZ GTO basis set ■ To approach the correct value of vertical 

266 attachment energy we changed the exponent of the uncontracted p-orbital of chlorine from 

267 0.162 to 0.094. In addition to this modification of the exponent, we also restricted the CAS 

268 CI model to four active electrons from the highest occupied molecular orbital (HOMO) and 

269 HOMO-1 and allowed their excitations into 11 lowest VOs. Using the notation of the Cg 

270 point group, this CAS model can be expressed as follows: 

(la' . . . 16a')''(la" . . . 7a")'^(17a' . . . 25a', 8a" . . . lla")^ (11) 

271 or 

(la' . . . 16a')3='(la" . . . r&y^Ua . . . 24a', 8a" . . . 12a")^ (12) 

272 as the ordering of VOs changes with nuclear geometry. Therefore, in this model we extend 

273 the space of VOs and reduce the number of active electrons as compared with model 1. 

274 Properties of the target calculated at the equilibrium nuclear geometry are summarized in 

275 Table [T] (we refer to this model as to Model 2 in the table and elsewhere in the text below). 

276 It shows that although the ground state energy is higher than the value calculated using 

277 model 1, the vertical excitation energy into the state is lower and closer to the published 

278 reference value. This CAS model of the target also gives lower vertical excitation energy to 

279 the state than model 1. Table [T]also shows that this modified model gives larger value 

280 of the dipole moment than model 1, however its value is still low enough to not to introduce 

281 a significant error in the position and width of the resonance. We also calculated the static 

282 dipole polarizability at equilibrium nuclear geometry using this model and found that its 

283 value is significantly smaller than that given by experiment. This can be due to the limited 

284 representation of the polarization effects by the primary GTO basis set. Although both our 

285 CAS models of the target have their limitations, model 2 using the modified primary GTO 

11 



286 basis represents correctly the vertical attachment energy, the quantity which is essential for 

287 the DEA calculations. 

288 In addition to the CAS model of the target we need a representation of the CF3 fragment 

289 in our LCP calculations in order to extrapolate the potential energy surfaces properly. In 

290 order to keep the model of the fragment consistent with the target model 2, we treated the 

291 CF3 radical at the SCF level as all the excitations of the active electrons in model 2 of the 

292 CF3CI target describe the C-Cl chemical bond. Therefore, any CI excitation model of the 

293 fragment would introduce a correlation which is not explicitly included in the target model. 

294 In order to check the quality of this representation we calculated the harmonic frequency 

295 of the umbrella mode. Its comparison with previously published value is given in Table [Tll 

296 It shows that our calculated harmonic frequency is slightly higher than the experimental 

297 value, but they are in a good agreement to confirm that our model is a suitable choice for 

298 representation of the CF3 radical. 



299 B. Scattering model 

300 The K-matrix calculations were performed using a sphere with radius rn=15ao. Cor- 

301 responding continuum basis set was represented by single-center uncontracted GTOs with 



302 exponents optimized by the program GTOBAS [28|. Partial waves up to / = 3 (9s, 7p, 7d, 



303 7f) were used. The deletion threshold in the orthogonalization procedure for the continuum 



304 orbitals 15| was set to 5thr = 9 x 10~^. This value was found by performing calculations 

305 in the static exchange approximation. It gives a stable representation of the scattering 

306 continuum and does not show any problems related to linear dependence of the continuum 

307 orbitals. 

308 In order to obtain the position and width of the fixed-nuclei resonance as a function of the 

309 nuclear geometry we performed a scattering calculation of the eigenphases and fitted them 

310 using the Breit-Wigner formula with the energy dependent background [l^ as discussed 

311 above. We tested both CAS CI models of the neutral target discussed above. First, the 

312 lowest 16 CI target states calculated using Model 1 were used in the CC expansion and 

313 the J?-matrix calculation at equilibrium nuclear geometry. This calculation gives converged 

314 results with respect to the number of target states included. 

315 The corresponding cross section for elastic "^Ai scattering is presented in Fig. [1] where we 
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FIG. 1. Cross section for the elastic electron scattering off CF3CI (^Ai symmetry) calculated in the 
fixed-nuclei approximation at equilibrium nuclear geometry. Calculations employing target model 



1 is compared with results calculated using model 2, with experimental results 

j 

R-matrix calculation ^291] . 



30| and with other 



318 also compare it with calculations of Beyer et al. |29| . This curve (designated as Target model 

318 1) shows the peak at energy around 3.3 eV that is significantly above the experimental value 

319 due to experimental work by Underwood-Lemons et al. 30| and the theoretical calculation 

320 of Beyer et al. [29]. Since a further increase of the number of target states included in 

321 the CC expansion did not lead to any considerable shift of the peak towards lower ener- 

322 gies, it suggests that this artificially high position of the resonance is not due to incorrect 

323 representation of the polarization effects j3l|, but rather due to an incompleteness of the 

324 target model. Therefore, we performed another set of scattering calculations using target 

325 CI model 2 described above. Since this model was constructed to reproduce the correct 
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326 vertical attachment energy, it also gives the correct position of the resonance peak in the 

327 R- matrix calculation carried at the equilibrium nuclear geometry (as plotted in Fig. [1]). The 

328 rapid increase of the cross section at lower energies is due to the long-range dipole potential 

329 included in the outer region. For a molecule with a nonzero permanent dipole moment and 



330 fixed orientation the elastic cross section diverges at zero energy [32|, |33| , the feature which 

331 is observed in our cross section, but not in the calculation of Beyer et al. 29|, apparently 

332 because the dipole effects were not included completely there. On the other hand, since 

333 our dipole moment at the equilibrium internuclear separation is too big (see Table I), it is 

334 evident that our elastic cross sections are strongly overestimated at low energies. 

335 There were 12 target states included in our CC expansion (three lowest in singlet and 

336 triplet state of A' and A" symmetries). The R- matrix was propagated in the dipole potential 

337 given by the target CAS CI model (as discussed above) and in the potential given by the 



338 dipole and quadrupole coupling of different scattering channels [15| . On the other hand, the 



339 R-matrix calculation of Beyer et al. [29| was performed at the level of static exchange with 

340 polarization which treats the electron correlation in a different way than our CAS CI model. 

341 This can be a partial reason of the quantitative difference between the two calculations at 

342 higher energies. 

343 Fig. [1] also presents the experimental measurements of the total elastic cross sections 



344 by Underwood-Lemons et al. [30|]. However, the problem with this and two other |12l . |34 | 

345 measurements for CF3CI is that it is not quite clear what is measured there. The total 

346 elastic cross section for a polar symmetric top is divergent even if rotations are included 



347 [35|]. Only the inversion splitting, which is extremely small for CF3CI, makes the elastic 

348 cross section finite [sG^. This means that the scattering amplitude at small scattering angle 

349 9 behaves like 1/9, and only at an extremely small angle 9inv it becomes finite. In the 

350 experiment of Underwood-Lemons et al. |3o| the elastic cross section is determined from 

351 the transmitted current under the assumption that transmitted electrons are not scattered, 

352 therefore what is measured in fact is the cross section integrated from a small angle 6'min 

353 to 180°, where 6'inin is determined by the geometry of the experimental apparatus, and 9m\r, 



354 is, most likely, significantly greater than 9inv In the experiment of Mann and hinder [12 1 

355 the total cross section is obtained by extrapolating the measured differential cross section 

356 to = 0. This procedure also gives an underestimated total cross section since the actual 

357 differential cross section behaves as 1/9'^ at small angles, if 6* > 9inv We think, therefore 
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358 that comparison between the theory and the experiment for the total (integrated) cross 

359 section is meaningless unless the angle Ojnm can be found from the experimental geometry. 

360 The purpose of plotting the experimental cross section in Fig. [T] is to demonstrate that our 

361 calculated resonance contribution is consistent with experimental results. 

362 At larger C-Cl separations, where the resonance turns into a bound state, the corre- 

363 sponding anionic bound state energies were calculated by diagonalizing the + 1 electron 

364 Hamiltonian constructed using the same target and scattering CAS CI model as was used 

365 for calculation of the resonance, just the integrals involving continuum GTOs were not re- 



366 stricted to the inner region of the sphere 



15|. This allows a good representation of the 



367 diffuse character of the anionic bound state with energy close to the autodetachment limit. 

368 The position and width of the fixed-nuclei resonance as a function of the nuclear geometry 

369 were obtained by fitting the eigenphases at energies close to the resonance to the Breit- 

370 Wigner formula with energy-dependent background performed at each nuclear geometry of 

371 our interest. The energy dependence of the background is predominately determined by 

372 the dipole moment of the neutral target. The corresponding value was obtained from the 

373 calculation of the target properties for every nuclear geometry and used as a parameter in 

374 the fitting procedure for this geometry. Fitting of the eigenphases calculated at equilibrium 

375 geometry using target model 2 gives a resonance position 1.954 eV, which is in a good 



correspondence with experimental value 1.83 eV due to Aflatooni and Burrow 



25|. 



377 IV. CALCULATIONS OF THE NUCLEAR DYNAMICS 

378 A. Dissociative electron attachment 

379 The total DEA cross section calculated using our two-dimensional LCP approach is plot- 
38H ted in Fig. [2l This graph shows its comparison with previously published one-dimensional 



382 non-local calculations 



13j. Our previously published two-dimensional LCP calculations 



383 give significantly larger magnitude of the total DEA cross section than experimental 



384 works 



25 



37( 1 and one-dimensional non-local semiempirical calculation 13|. This prob 



386 lem was attributed to incorrect dependence of the width function (which was arbitrarily 

386 extended to two dimensions) on the coordinate r. As can be seen in Fig. [21 results of our 

387 present two-dimensional LCP calculations using the complex PES constructed from the 
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FIG. 2. Total DEA cross section calculated using the model discussed in the present work is 
compared with experimental results of Underwood-Lemons et al. [37], Aflatooni and Burrow {25! 
and with previously published one-dimensional non-local calculations due to Wilde et al. [3]. 



3B8 J?-matrix results are in very good agreement with t 



le non-local calculations as well as with 



389 experimental results due to Aflatooni and Burrow 25|]. Although our total cross section is 

390 smaller than the experimental results of Underwood-Lemons et al. [sT^, it is closer to the 



391 measurements of Aflatooni and Burrow 



25|, and the position of the peak is in very good 



392 correspondence with this experimental work. All this suggests that the two-dimensional 

393 PES constructed from the results of the R-matrix calculations has more correct dependence 

394 on both reaction coordinates than the previously published model [li]. It results in a correct 

395 shift of the peak in the total DEA cross section with respect to the VAE as well as in the 
398 correct magnitude of the cross section. One additional note regarding the complex PES 

398 should be made here. The DEA cross section is related to the position of the crossing seam 

399 between the neutral and anionic PES. Its shape and position with respect to the minimum 
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FIG. 3. The crossing seam between the neutral and anionic two-dimensional PES obtained from 

n 

the present Ji-matrix calculations and its comparison with model potential used in pj]. The cross 
(circle) denotes the equilibrium geometry of the neutral PES obtained in the present work (used 
in Q). 



400 of the neutral PES determine where the anionic system becomes stable. Subsequently, it 

401 has an influence on the survival probability. The crossing seams obtained from the present 

402 J?-matrix calculations and from the model potential [l| are compared in Fig. |3l The relative 

403 position of the crossing seam is similar in both models. This supports our argument that 

404 the difference in magnitude of the DEA cross section between our present model and the 

405 results published in [l| is mainly due to the different width function F(p, r), rather than due 

406 to the substantial difference in the crossing seam between these two models. The oscillatory 

407 structure in the crossing seam obtained from the present R-matrix calculations is mainly 

408 an artifact of too low density of the grid of nuclear geometries used to calculate the PES. 
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409 In addition, the fact that the complex energies in the region, where the anion is metastable 

410 are obtained in different way than the eneries of the bound anion, also raises the numerical 

411 issues with exact determination of the crossing seam and partially also contibutes to the 

412 oscillatory structure present in Fig. [31 



413 Comparisons of local two-dimensional and one- dimensional results with nonlocal one- 

414 dimensional results are given in Refs. The agreement between the local and non- 
415 local results is very good because the resonance occurs at a relatively large energy, and the 
416 long-range interaction plays a minor role. 



417 Fig. m shows the distribution of different final vibrational states of the CF3 fragment 

418 calculated using the complex PES constructed from the R- matrix results. In Fig. [5] we show 

419 the distribution calculated using the PES obtained in Ref. l|] but corrected in such a way 

420 that PES in the intermediate region smoothly turns into the PES in the asymptotic region 

421 avoiding the mismatch discussed in Sec. [TTl Both graphs show that the low vibrationally 

422 excited states of the CF3 fragment will be more populated than the ground state and in 

423 both figures the positions of the peaks rise with increasing quantum number u. However, 

424 each calculation predicts highest population for a different excited state. While our complex 

425 PES constructed from the R- matrix results gives the highest peak for u = 1, the PES 

426 described in ij leads to the highest peak for vibrational state z/ = 3, as can be seen in 

427 Fig. O Total DEA cross section is mainly determined by the complex PES in the region 

428 where the anionic system is not bound and does not strongly depend on the behavior of 

429 the PES in the region where the negative ion is stable. On the other hand, the distribution 

430 of vibrational states of the fragment is strongly influenced by the final-state interaction in 

431 the region where the anion is stable. In both calculations shown in Fig. H] and Fig. [5], this 

432 region of the PES was partially modeled (as discussed above and in to achieve correct 

433 asymptotic behavior of the potential. Although our PES constructed from the R-matrix 

434 results is free of several limitations of the surface described in the extrapolation of the 

435 bound-state anionic surface in both cases makes it difficult to decide, how quantitatively 
438 reliable this distribution is in both calculations. 



18 



I 

o 




FIG. 4. DEA cross sections for different final vibrational states of the fragment CF3 calculated 
using the complex PES constructed from our ii-matrix results. 

B. Vibrational excitation 



440 Results of the resonant VE calculations using the PES based on the R-matrix results 

441 are plotted in Fig. |6l This figure shows the cross section for VE of the target from the 

442 ground state to the lowest excited state of the C-Cl stretching mode (denoted as (0, 1) in 

443 the figure) and the lowest excited state of the C-F deformation mode (denoted as (1,0) in 

444 the figure). These graphs show a very good agreement with our previous calculations [l| 

445 for both vibrational modes. However, our calculations agree with experimental results by 

446 Mann and hinder 12| for the umbrella mode, while both our models lead to approximately 



447 three times smaller magnitude of the cross section for the C-Cl stretching mode excitation 

449 than that measured by Mann and hinder I12I . 

n 

450 The previously published one- dimensional non-local calculation of VE [13\ takes into 
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FIG. 5. DEA cross sections for different final vibrational states of the fragment CF3 calculated using 
the model complex PES published previously [l| after the correction of the asymptotic behavior of 
channel potentials. 



451 account the C-Cl stretching vibrational mode only and the width function T{R) was adjusted 

452 to give cross sections corresponding to experimental results for this stretching vibrational 

453 modes excitation due to Mann and Linder [l^ . The low values of the present cross sections 

454 can be partially explained by introduction of additional channels of vibrational excitation 
456 in our two-dimensional calculations and lower flux towards the channel (0, 1). 

456 Although Fig. E] shows that the model constructed using the R- matrix results leads to VE 

457 cross sections very similar to those calculated in Ref. jj, Fig. [7] shows that these two models 
45B predict different results for higher final vibrational states. To our best knowledge, there 

459 are no experimental data for VE to these states to compare with our cross sections. The 

460 fact that the differences between our two models become more significant with increasing 
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E (eV) E (eV) 

FIG. 6. VE cross sections from the ground vibrational state to the lowest excited state of the C-Cl 
stretching mode (a) and the lowest excited state of the umbrella mode (b). Results calculated 
using the R-matiix complex PES are compared with cross sections obtained from the previously 

n n 

published model and experimental results due to Mann and Linder [1^ . 

461 final vibrational state is understandable, since the target vibrational eigenf unctions become 

462 spatially more extended with increasing vibrational state {1^2, 1^3) ■ This means that the 

463 different behavior of corresponding complex PES farther from equilibrium will have larger 
468 influence on results, as can be seen in Eq. @. 



466 V. CONCLUSION 

467 The ab initio R-matrix method allowed us to calculate the complex two-dimensional 

468 PES for CF3CI collisions. We used then DVR method to obtain the solution of coupled 
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FIG. 7. VE cross sections for excitation from the target vibrational ground state to higher excited 
states. Results obtained from the complex PES constructed from the fi-matrix results (solid line) 



are compared with cross sections calculated using the PES described in 



1[ (dashed line). 



469 stationary equations dlj, and DEA and VE cross sections for the e— CF3CI collision process. 

470 Our results for the total DEA cross section and the cross section for VE of the umbrella 

471 mode agree quite well with experiments, In addition we obtained the final-state vibrational 

472 distribution in the CF3 fragment free of instabilities found in our previous calculations. 

473 However our cross sections for VE of the C-Cl stretching mode are significantly lower than 

474 the experimental results of Mann and Linder [l2| and the results of previous semiempirical 

475 one-dimensional calculations This is something one might expect because of the extra 

476 inelastic channels of excitation of umbrella mode that leads to a redistribution of the flux. 

477 However, disagreement of the present 2-dimensional results with the experiment is puzzling 
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478 and requires a further investigation. 
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